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ABSTRACT 

We  present  recent  developments  in  the  theory  of  a 
new  class  of  active  multivariate  FT  spectrometers  and 
describe  their  potential  uses  in  the  optical  detection  and 
identification  of  improvised  explosive  devices  (IEDs)  and 
CB  agents  and  contaminants.  This  technology  is  based  on 
the  use  of  dual  ultra-high-speed  (UHS)  birefringence 
interferometers  for  the  simultaneous  and  independent 
modulation  of  both  illumination  and  collected  beams,  and 
employs  photoelastic  modulators  (PEM)  as  the 
birefringent  elements  in  these  interferometers.  With 
typical  scan  rates  of  104-105  interferograms  per  second 
and  broad  spectral  range,  these  instruments  can  rapidly 
collect  multivariate  spectral  data  describing,  in  a 
correlated  way,  the  comprehensive  excitation,  emission, 
and  delay/lifetime  spectral  properties  of  the  sample/target. 
The  intrinsic  process  of  data  fusion  that  melds  the 
multivariate  spectral  data  into  one  interferogram  signal 
allows  for  the  simultaneous  and  correlated  use  of  all  three 
spectral  variables  for  improved  discrimination  and 
identification  of  chemical  species  of  interest. 

1.  INTRODUCTION 

PEM-based  FT  spectrometers  can  easily  acquire  105 
spectra  (interferograms)  per  second  and  cover  broad 
spectral  ranges  (UV  to  NIR  with  fused  silica  optics,  and 
NIR  to  mid-IR  with  ZnSe  optics).  PEM-based 
interferometers  can  be  used  simultaneously  to  modulate 
both  the  illumination  beam  and  the  collected  light.  In  such 
a  system,  the  chirped  modulation  (Stann  et  al.,  1996)  of 
the  illumination  beam  (Buican,  1993,  1995)  encodes  a 
time  origin  for  each  half-period  of  the  PEM,  as  well  as  the 
spectrum  of  the  illumination  beam.  The  second  PEM- 
based  interferometer  additionally  encodes  the  spectrum  of 
the  collected  light.  The  resulting  interferogram  represents 
a  rich  fusion  of  correlated  data  and  presents  the  promise 
of  improved  detection,  discrimination,  and  identification 
of  chemical  species  through  digital  postprocessing. 

We  present  in  this  paper  a  powerful  and  elegant 
mathematical  interpretation  of  the  linear  mapping  of 
multivariate  PEM/FT  spectrometry  (M-PEMS).  This 
interpretation  exposes  the  way  in  which  M-PEMS  is  a 


direct  generalization  of  standard  FT  spectrometry  by 
separating  the  complex  effects  sample  (specific 
multivariate  spectrum)  and  instrument  (operating 
parameters  and  scanning  characteristics)  have  on  the 
interferogram.  Furthermore,  it  sheds  light  on  how  an  M- 
PEMS  instrument  should  be  designed  and  operated. 

We  further  show  that  the  separation  into  sample  and 
instrument  components  leads  to  an  understanding  of  M- 
PEMS  systems  as  data  fusion  optical  processors,  and  of 
how  this  important  function  can  be  externally  controlled. 
This  further  opens  the  way  to  adaptive  operation  for 
improved  spectral  discrimination  and  identification  of 
certain  threats. 

Finally,  we  discuss  the  strategic  advantages  offered 
by  multivariate  PEM/FT  regarding  (i)  the  protection  of  U. 
S.  and  Allied  troops  against  IEDs  and  (ii)  the  early 
warning  of  dangerous  environmental  CB  threats  evolving 
in  the  open  battlefield  or  in  the  homeland.  The  design  and 
development  of  remote  and  standoff  sensing  systems 
based  on  M-PEMS  may  include  lidars,  ladars,  and 
hyperspectral  imagers.  All  these  systems  will  greatly 
benefit  from  the  massive  data  fusion  intrinsically 
performed  by  the  optical  systems  of  M-PEMS  sensors. 
Such  sensors  may  thus  even  provide  the  basis  for  very  fast 
automated  defense  system. 

2.  PEM/FT  SPECTROMETRY  FUNDAMENTALS 

2.1  UHS  Spectral  Modulation  Using  PEM-Based 
Interferometers 

PEM/FT  spectrometers  are  based  on  birefringence 
interferometers,  which  introduce  a  variable  retardation 
between  interfering  waves  separated  by  their  mutually 
orthogonal  directions  of  polarization.  The  retardation  is 
proportional  to  the  instantaneous  birefringence  of  the 
modulator  and  can  be  varied  in  time  in  order  to  scan  an 
interferogram.  PEM-based  FT  spectrometers  cover  broad 
spectral  ranges  (UV  to  NIR  with  fused  silica  optics,  and 
NIR  to  mid-IR  with  ZnSe  optics).  In  its  simplest 
implementation  (Fig.  1),  a  birefringence  interferometer 
consists  of  a  birefringent  modulator  placed  between  two 
polarizers  (Buican,  1985,  1990).  The  optical  axes  of  the 
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birefringent  modulator  are  at  45°  relative  to  the  direction 
of  the  polarizers. 


Fig.  1:  The  basic  birefringence  interferometer.  E^,  Eout,  input 
and  output  electric  vectors;  M,  birefringence  modulator;  nx,  ny, 
principal  indices  of  refraction  of  the  birefringence  modulator 
and  the  corresponding  optical  axes;  and  P1?  P2,  polarizers  and 
their  optical  axes. 

The  birefringence  interferometer  in  Fig.  1  can  be 
regarded  as  an  optical  frequency-dependent  intensity 
modulator,  with  a  spectral  modulation  function 
MAn(a))  =  lout(o))/lin(o )  controlled  by  the  instantaneous 

birefringence  A n=nx-ny.  The  amplitude  of  the  phase 
difference  for  interfering  monochromatic  waves  of 
frequency  a>  is  ®  =  yco l,  where  the  retardation  amplitude 
y-  Tnax  is  the  maximum  time  delay  between  the  interfering 
waves.  For  a  lossless  interferometer,  the  modulation 
function  is  mt(©)  =  [i  +  cos(t©)]/2  .  For  a  birefringence 
interferometer  using  PEMs,  the  temporal  retardation  is 
given  by  x(t)  =  ycosin(co0t) ,  where  co0  is  the  angular 
frequency  of  the  PEM’s  mechanical  oscillation,  and  the 
retardation  amplitude  is  given  by  y  =  Bl/c,  where  B  is  the 
birefringence  amplitude,  and  /  is  the  thickness  of  the 
birefringent  modulator.  Thus,  for  PEM-based  birefringent 
interferometers,  the  spectral  modulation  function  is 

M  (co,  t )  =  (co)  /  Iin  (co)  =  [l  +  b  (yco,  ft)0r)] ,  ( 1 ) 

where  b(®,y/)  =  cos(®  siny/) .  The  harmonic  rather  than  linear 
time  dependence  of  the  interferogram  phase  for  a  PEM- 
based  interferometer  at  first  appears  to  be  an 
inconvenience  but  in  fact  plays  an  essential  role  in  active 
multivariate  spectrometry. 

2.2  Active  PEM/FT  Spectrometry 

We  developed  the  concept  of  multivariate  PEM/FT 
spectrometry  (Buican,  1993,  1995),  in  which  both  the 
illumination  beam  and  the  collected  light  pass  through 
PEM-based  interferometers  (Fig.  2).  In  such  a  system,  the 
chirped  modulation  of  the  illumination  beam  encodes  a 
time  origin  for  each  half-period  of  the  PEM,  while  also 
encoding  the  spectrum  of  the  illumination  beam.  This 
information  can  subsequently  be  used  in  order  to  extract 


from  the  total  collected  intensity  waveform 
(interferogram)  the  delays  and  delay  distributions 
associated  with  each  illumination  wavelength.  By  also 
passing  the  collected  light  through  a  PEM-based 
interferometer,  one  can  further  encode  into  the 
interferogram  the  spectrum  of  the  collected  light.  The 
latter  can  differ  from  the  illumination  spectrum  if  inelastic 
scattering  and/or  fluorescence  are  present. 


Fig.  2:  Conceptual  design  of  a  trivariate  PEM/FT  spectrometer 
operating  in  the  infrared.  A,  input  aperture;  B,  beam  splitter;  CS, 
wavelength  calibration  source;  CD,  wavelength  calibration 
detector;  D,  detector;  DM,  demultiplexer;  M,  off-center 
parabolic  mirror;  Pin,  Pout,  input  and  output  polarizers;  and 
PEMs,  photoelastic  modulators.  The  illumination/excitation 
beam  is  represented  in  yellow,  while  the  collected/emitted  light 
is  represented  in  red. 

The  interferogram  generated  by  the  detector  of  a 
multivariate  spectrometer  is  a  linear  transform  of  an 
intensity  spectrum  which  itself  is  a  function  of  three 
independent  variables:  the  illumination  frequency,  the 
collection  frequency,  and  the  delay  between  illumination 
and  collection.  In  a  fluorescence  microscope  or  flow 
cytometer  (Fig.  3a),  the  only  significant  delays  between 
the  modulation  of  the  illumination  beam  and  that  of  the 
collected  light  are  directly  related  to  fluorescence  lifetime 
and  the  interferogram  generated  by  the  dual  PEM/FT 
spectrometer  encodes  a  sampling  of  the  FT  of  an  intensity 
spectrum  that  that  is  simultaneously  a  function  of 
excitation  optical  frequency,  emission  optical  frequency, 
and  fluorescence  lifetime.  In  a  lidar  application  (Fig.  3b), 
the  propagation  delays  dominate  and  the  multivariate 
spectrum  involved  is  a  function  of  excitation  and  emission 
optical  frequencies,  and  of  propagation  time  (range). 

3.  THEORY  OF  MULTIVARIATE 
SPECTROMETRY 

3.1  The  Direct  Problem 

The  interferogram  produced  by  a  multivariate  system 
such  as  the  lidar  in  Fig.  3b  is  a  linear  transformation  of  a 
“trivariate  spectrum”,  which  is  the  distribution  of 
collected2  light  intensity  simultaneously  over  excitation 
frequency,  emission  frequency,  and  delay  (range);  we 


2  The  effect  of  the  illumination  and  collection  geometries  and  optics  can 
be  factored  out  to  leave  a  “true”  trivariate  spectrum  describing  the 
intrinsic  spatial  and  optical  frequency  properties  of  the  target. 


1  In  the  context  of  multivariate  spectrometry,  we  choose  the  temporal 
rather  than  spatial  versions  of  frequency  and  retardation  because  one  of 
the  variables  of  interest  is  the  lifetime  or  propagation  delay. 
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regard  these  as  the  natural  spectral  variables.  If 
s(cox,com,p)  is  the  tri variate  spectrum  in  the  natural 


Fig.  3:  Multivariate  PEM/FT  spectrometry.  Simplified  diagrams 
of  two  previously  proposed  instrumental  configurations  for 
multivariate  spectrometry:  a  scanned  multivariate  hyper  spectral 
fluorescence  imager/microscope  (a);  and  a  multivariate 
fluorescence  lidar  (b).  DM,  dichroic  mirror3;  Pin,  Pout,  input  and 
output  polarizers;  PEM,  photoelastic  modulator;  L,  lens 
(microscope  objective  optics). 

variables,  then  the  interferogram  is  given  by  the  following 
integral  transform: 

/„(#)  =  «J J J  [i+b (r,fflx,ffl0,(;-p)+p0)]x 

0  0  0  ’ 
x  [l  +  b  {ymcom ,  cQqJ)  ]S'(a>x,  com ,  p) dcoxdcomdp 

where  S'(cox,com,p)  =  Is(cox)S(cox,a>m,p)R(com)  is  the 
effective  tri  variate  spectrum,  with  Is(co)  the  intensity 
spectrum  of  the  light  source  and  R(co)  the  spectral 
responsivity  of  the  detector;  yxm  are,  respectively,  the 

(temporal)  retardation  amplitudes  of  the  excitation  and 
emission  interferometers;  co0x,co0m  are,  respectively,  the 
excitation  and  emission  interferometers’  angular 
frequencies  of  oscillation;  <p0  is  the  phase  difference 
between  the  excitation  and  emission  interferometers;  and 
a  is  a  system-dependent  coefficient.  This  is  the  direct 
problem  of  trivariate  PEM/FT  spectrometry. 

It  is  obvious  from  Eq.  (2)  that  the  transform  mapping 
the  trivariate  spectrum  to  the  interferogram  is  related  to 
the  cosine  Fourier  transform  and  will  thus  impose  certain 
symmetry  constraints  on  the  system’s  spectral  response.  It 
is  also  obvious  that  it  maps  functions  of  three  variables 

3  A  dichroic  mirror  is  used  in  fluorescence  applications  to  separate  the 
emission  and  excitation  wavelengths.  For  lidar  applications  in  which 
backscattered  light  is  important,  a  duplexer  mirror  will  be  used  instead. 


(the  trivariate  spectrum  S'(cox,com,p))  to  functions  of  one 
variable  (the  interferogram)  and  therefore  cannot  be 
inverted  in  the  strict  sense  of  the  word.  Additionally,  the 
retardation  in  PEM-based  interferometers  is  a  harmonic 
function  of  time.  Finally,  a  dual  interferometer  system  has 
five  adjustable  parameters4  that  affect  the  interferogram. 
Together,  these  factors  seem  to  make  multivariate 
PEM/FT  interferograms  nearly  intractable.  However,  we 
show  in  the  next  subsection  that  there  is  an  elegant 
geometric  interpretation  of  the  direct  transform  whereby 
the  interferogram  can  be  obtained  by  sampling,  along  a 
curved  3-D  path,  a  function  of  three  variables  which  itself 
is  related  to  S'  (cox ,  com ,  p)  through  a  linear  transform. 

3.2  Separation  of  Sample  and  Instrument 
Descriptions 

As  the  transform  in  Eq.  2  is  linear,  we  will  initially 
concern  ourselves  with  a  point  spectrum  in  the  space  of 
natural  coordinates,  located  at  excitation  and  emission 
frequencies  ux  ,  um  ,  and  delay  v,  and  which  describes  a 
target  that  is  excited  only  at  an  optical  frequency  ux  and 
which  subsequently  emits  at  the  frequency  um  after  a 
delay  v, 

S'  (cox,com,  p)  =  S(cox  -ux)S(com  -um)S(p-v)  (3) 

(Fig.  4a).  From  Eq.  (2),  the  interferogram  for  this  point 
spectrum  is 

ld(t)  =  a[\+b ( yxux , co0x (t - v) +  % ) ][l  +  b  ( ymum , co0mt) ] .  (4) 

If  we  change  from  the  natural  coordinates  {cox,com,p}  to 
the  homogeneous  ones  {cox,co2,co 3}  in  which  the  natural 
coordinates  become  cylindrical  (Fig.  4b)  by  converting 
the  delay  coordinate  p  to  a  phase  delay  angle  relative  to 
the  phase  of  the  excitation  interferometer, 

cox  =  cox  cos (&>0xp),  co2  =  cox  sin(&>0xp),  co3  =  com  ,  (5) 

the  original  effective  tri  variate  spectrum  S'(a>x,com,p ) 
yields  a  new  tri  variate  spectrum  Sh'{cox,co2,a>3)  in  the 
new,  homogeneous,  variables.  We  further  construct  a 
symmetrized  spectrum  S'hs  by  adding  to  S'h  its 
reflections  in  the  { colyco2}  plane  and  the  co3  axis5  and 
finally  we  augment  it  to  S' hsa  by  adding  projections  onto 
the  {a>x,co2}  plane,  the  co3  axis,  and  the  origin6  (Fig.  4b). 

A  3D  Fourier  transform  ^  of  the  new  spectrum 
S' hsa  yields  a  tri  variate  “interferogram”  I  d3  =  ^S’hsa , 
which,  because  of  the  symmetry  properties  of  this 

4  Two  retardation  amplitudes;  two  scanning  frequencies;  and  one  relative 
phase. 

5  The  symmetries  added  to  the  trivariate  spectrum  are  a  generalization  of 
the  inversion  symmetry  in  standard  FT  spectrometry. 

6  The  projections  added  to  the  tri  variate  spectrum  correspond  to  the 
effects  of  the  d.c.  terms  in  the  interferogram  of  Eq.  2. 
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spectrum,  is  a  real  function  of  three  time -like  variables 
(temporal  retardations),  ldt3(r)  ,  where  t  =  is  a  3D 


Fig.  4:  Tri variate  point  spectrum  in  natural  coordinates  (a)  and 
homogeneous  frequency  coordinates  (b).  The  point  spectrum  in 
(b)  has  been  symmetrized  and  augmented.  um,  ux,  v  are, 
respectively,  the  emission  and  excitation  optical  frequencies, 
and  the  delay  between  excitation  and  emission  for  the  trivariate 
point  spectrum;  w,  unit  of  weight  for  the  point  distributions  in 
space  of  homogeneous  frequency  coordinates;  \) /,  angular 
coordinate;  G)pl;_5,  position  vectors  in  the  space  of  homogeneous 
frequency  coordinates  for  the  various  point  distributions;  co0x, 
mechanical  oscillation  frequency  of  the  excitation  PEM-based 
interferometer. 


“time”  vector.  For  the  symmetrized  and  augmented  point 
spectrum  in  Fig.  4b,  the  resulting  3D  interferogram  is 

Id,  3  W  =  y5[2cos(®^r)+  2cos(a/4r)+  2cos(®^r)+ 

+  c°s(®^r)+c°s(fi/2r)] 


where  p  is  a  constant  coefficient  and  the  frequency 
vectors  correspond  to  the  positions  shown  in  Fig.  4b.  This 
is  the  “trivariate”  interferogram,  a  real  function  over  the 
space  of  3D  time  vectors  conjugate  to  the  frequency 
vectors  defined  in  Eq.  (5).  It  is  the  sampling  of  this 
trivariate  interferogram  along  a  path  parametrized  by  the 
actual  time,  p  =  {r(t)},  that  results  in  the  interferogram 
generated  by  the  detector,  Id(t)  =  Idt3(r(t))  •  For  the  system 

that  generates  the  interferogram  given  by  Eq.  (4),  the 
sampling  path  can  be  shown  to  be 


r,0) 

rxsm{co0xt  +  %x) 

r2(0 

= 

-yxcos{co0xt  +  %x) 

t3(0_ 

7,„  sin(®0„,0 

As  discussed  in  subsection  3.3,  this  is  a  Lissajous  figure 
on  a  cylindrical  surface.  The  time-dependent  intensity 
resulting  from  the  sampling  of  the  trivariate  interferogram 
in  Eq.  (6)  along  the  path  in  Eq.  (7)  can  be  obtained  by 
substituting  in  the  former  the  time  vector  from  the  latter. 
The  result  can  be  shown  to  be  the  expected  ID  time- 
dependent  interferogram  given  by  Eq.  (4)  for  the  trivariate 
point  spectrum  in  Eq.  (3). 

It  is  important  to  notice  at  this  point  that  the  3D 
interferogram,  exemplified  for  the  case  of  a  point 
tri  variate  spectrum  by  Eq.  (6),  is  independent  of 


instrument  parameters  and  exclusively  describes  the 
tri  variate  spectral  properties  of  the  sample/target.  Because 
of  the  linear  nature  of  the  system,  the  same  will  be  true 
for  any  arbitrary  superposition  of  point  spectra,  and,  thus, 
for  arbitrary  trivariate  spectra.  On  the  other  hand,  the 
sampling  path  given  by  Eq.  (7)  only  depends  on  the 
instrument’s  parameters,  such  as  the  retardation 
amplitudes  of  the  two  PEM-based  interferometers,  their 
mechanical  oscillation  frequency,  and  their  relative  phase. 
Then,  the  actual  ID  detector  interferogram  for  any 
trivariate  spectrum  can  be  obtained  by  the  same  process 
of  sampling  the  corresponding  trivariate  interferogram 
along  the  instrument-determined  sampling  path  of  Eq.  (7). 

Thus,  while  the  multivariate  interferogram  is  a 
description  of  the  complex  spectral  properties  of  the 
sample,  the  geometry  of  the  curved  path  is  exclusively 
determined  by  the  instrument's  operating  parameters  and 
scanning  characteristics.  The  two  domains  are  reunited  by 
sampling  the  multivariate  interferogram  along  the 
instrumental  path,  which  results  in  the  actual  time- 
dependent  interferogram  as  generated  by  the  detector. 
Resolving  the  complex  trivariate  spectrometry  problem 
into  a  straightforward  3D  FT  and  the  geometry  of  the 
sampling  path  conveniently  separates,  conceptually  and 
computationally,  the  spectral  properties  of  the  sample 
from  the  properties  of  the  instrument. 

This  interpretation  also  exposes  the  way  in  which 
multivariate  PEM/FT  spectrometry  is  a  direct 
generalization  of  standard  FT  spectrometry.  Indeed,  in 
standard  FT  spectrometry,  the  detector  signal  is  a  function 
of  retardation,  from  which  the  actual  interferogram  signal 
is  derived  by  scanning  the  retardation  as  a  function  of 
time.  In  the  multivariate  case,  the  detector  signal  is  still  a 
function  of  retardation,  which  is  now  a  3D  vector,  and  the 
detector  signal  is  again  obtained  by  scanning  the  3D 
retardation  as  a  function  of  actual  time,  i.  e.,  along  a  path 
in  3D  space. 

By  applying  an  inverse  3D  FT  transform  to  both  the 
trivariate  interferogram  and  the  sampling  path  results, 
respectively,  in  the  (symmetrized  and  augmented) 
tri  variate  spectrum  S' hsa  on  the  one  hand,  and  the  system’s 
point  response  function  on  the  other.  This  is  the  3D 
frequency  space  equivalent  of  the  separation  of  the 
problem  into  sample  and  instrument  specific  parts,  with  a 
convolution  of  the  spectrum  and  point  response  function 
as  the  equivalent  of  the  sampling  along  a  path  in  the  space 
of  time-like  coordinates. 

3.3  Geometry  of  Interferogram  Sampling 

The  trivariate  interferogram,  which  is  a  function  of 
three  time-like  variables,  cannot  be  arbitrarily  sampled. 
However,  the  temporal  evolution  of  the  phases  of  the  two 
interferometers  results  in  a  ID  sampling  path  in  the  space 
of  these  three  time-like  variables.  When  sampled  in  this 
way,  the  trivariate  interferogram  yields  information  that 
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can  be  used  to  reconstruct,  in  a  unitary  and  homogeneous 
way,  the  original  tri variate  spectrum7. 

The  sampling  path  given  by  Eq.  7  is  a  Lissajous 
figure  (Bogle  et  al,  1994)  on  the  surface  of  a  cylinder 
centered  at  the  origin;  with  its  axis  parallel  to  the  r3  axis 
and  spanning  the  interval  t3  e  [-  ym,ym\ ;  and  of  radius  yx 
(Fig.  5).  We  call  this  the  sampling  cylinder.  The  origin  of 


Fig.  5:  Sampling  path  for  trivariate  interferograms.  (See  the  text 
for  definitions  of  the  various  parameters.) 

the  angular  coordinate  in  the  {tvt2}  plane  is  at  -n  12 ,  and 
the  origin  of  the  path  is  rotated  by  cp0x  relative  to  it.  The 
sampling  path  is  a  closed  curve  if  the  ratio  of  the  PEM- 
based  interferometer  frequencies,  ojQxlco0m ,  is  a  rational 
number.  In  the  simplest  case,  the  two  frequencies  are 
equal  and  the  path  is  an  ellipse  with  its  plane  rotated  by 
the  initial  phase  cp0x .  Several  closed  paths  with  zero  initial 

phases  are  shown  in  the  first  column  of  Fig.  6.  It  is 
obvious  that,  by  manipulating  the  ratio  of  the  two  PEM 
frequencies,  one  can  achieve  as  fine  a  scanning  of  the 
sampling  cylinder  as  desired,  but  at  the  cost  of  increased 
scanning  time.  Indeed,  if  co0x/co0m  =n/m ,  where  m  and  n 
are  coprime,  then  m  gives  the  number  of  full  oscillations 
of  the  sampling  path  along  the  r3  axis,  while  n  is  the 
number  of  full  circles  in  the  {r1?r2}  plane,  for  one  full 
period  of  the  closed  sampling  path.  Thus,  the  time 
required  to  travel  the  full  sampling  path,  and  thus  generate 
a  full  interferogram  and  acquire  all  the  accessible  spectral 
information,  is  Tscan  =2mnlco0m  and  increases  linearly  with 
the  density  of  the  coverage  of  the  sampling  cylinder  by 
the  sampling  path. 

In  the  limit  of  infinite  scanning  time,  the  sampling 
path  passes  through  every  point  of  the  sampling  cylinder. 
This  ideal  case  represents  the  best  one  can  do  in  terms  of 
sampling  the  trivariate  interferogram  for  given  retardation 
amplitudes  yx,ym .  The  ideal  spectral  point  response 
function  can  be  easily  derived  as  the  inverse  3D  FT  of  the 
sampling  cylinder,  and  is  given  by 

Hr„y.A®)  =  Kh(rxy[a>[+®~2  ) sinc^,,,^)  ,  (8) 

7  Within  the  resolution  allowed  by  the  spectral  response  function 
corresponding  to  the  instrument’s  operating  parameters. 


where  k  is  a  constant.  This  ideal  spectral  point  response 
function  is  illustrated  in  the  last  row  of  Fig.  6,  where  one 
can  notice  that  the  response  spreads  out  more  in  the 
excitation/delay  plane  than  along  the  emission 

axis  co3 ,  as  the  former  is  governed  by  the  Bessel  function 
J0,  while  the  later  by  a  sine  function.  One  can  also  see 
that,  for  n^m,  the  sampling  path  spectral  response 
functions  quickly  approach  that  for  the  sampling  cylinder 
itself.  Indeed,  even  for  n/m=6/ 5,  the  spectral  response 
function  in  a  plane  containing  the  emission  frequency  axis 
differs  from  the  ideal  one  mostly  in  the  higher-frequency 
diagonal  regions  of  the  plots  which  are  associated  with 
the  angular  periodicity  of  the  path  on  the  cylinder.  At 
n/m= 21/20,  the  angular  period  of  the  path  has  already 
decreased  to  the  point  where  these  diagonal  contributions 
have  moved  to  even  higher  frequencies  and  out  of  the 
plots. 

4.  DATA  FUSION  AND  POSTPROCESSING 
4.1  Intrinsic  Data  Fusion 

Data  fusion — the  process  of  combining  and 
representing  disparate  data  in  a  correlated  and 
homogeneous  form  that  can  be  effectively  and  efficiently 
processed  by  discrimination  and  identification 
algorithms — is  a  major  requirement  for  any  realistic  CB 
or  IED  warning  system.  It  allows  small  but  specific 
changes  in  individual  detector  signals  to  be  extracted  from 
the  noise  of  the  environment  by  identifying  specific 
correlations  that  are  absent  from  the  stronger  but 
nonspecific  noise.  In  the  M-PEMS  context,  we  refer  to 
this  as  the  Correlated  Multivariate  Advantage. 

Conventional  data  fusion  is  performed  numerically 
on  the  signals  collected  from  a  multitude  of  separate 
sensors  and  may  require  sophisticated  software  and 
considerable  data  processing  and  storage  resources.  M- 
PEMS  fuses  data  for  three  important  optical  spectral 
variables,  namely  excitation  and  emission  frequency  and 
emission  delay/lifetime  by  means  of  the  massively 
parallel  processing  power  of  its  optical  system,  and 
generates  one  signal,  the  interferogram,  that  encodes  the 
fused  data.  Thus,  the  task  of  digital  postprocessing  is  now 
reduced  to  discrimination  and  identification.  While 
nontrivial,  this  task  can  be  relatively  easily  performed  in 
compact,  high-performance  hardware  that  can  be 
embedded  into  a  mobile  and  even  portable  detection 
system. 

A  multivariate  PEM/FT  spectrometer  functions  as  a 
massively  parallel  optical  computer  which  computes  in 
parallel  at  all  excitation  and  emission  frequencies,  and 
time  delays.  Of  course,  a  standard  FT  spectrometer  also 
optically  computes  in  real-time  and  in  parallel  the  (cosine) 
FT  of  a  spectrum,  but  a  multivariate  PEM/FT  takes  this 
massively  parallel  computation  to  a  new  level  by  carrying 
out  at  full  speed  a  3D  FT  and  a  ID  sampling. 


The  three  types  of  data  optically  processed  and  fused 
by  an  M-PEMS  are  brought  together  in  the  homogeneous 


and  correlated  representation  of  the  trivariate  spectrum 
and,  through  an  optically  computed  3D  FT,  in 


n/m  Sampling  manifold  { 03b  co3 }  {oo2,  oo3}  { 03b  co2 } 


Fig.  6:  Trivariate  spectral  response  functions  for  various  sampling  paths.  The  first  column  shows  various  sampling  manifolds,  including 
both  sampling  paths  ( co0x  / co0m  =n/ m  values  of  1/1,  6/5,  and  21/20)  and,  in  the  limit,  the  full  sampling  cylinder.  The  other  three  columns 

show  three  orthogonal  sections  (in  the  planes  {co1,cd3},{co2,co3},{co1,co2}  )  through  the  corresponding  response  functions.  The  horizontal 
and  vertical  axes  of  the  response  function  plots  correspond  respectively  to  excitation/delay  and  emission  frequencies  ( cox  =  (cof  +  co\ )1/2 , 
com  =co3.).  These  plots  extend  to  approximately  10  resolution  widths  from  the  origin,  which  is  at  the  center  of  each  plot.  Similar  results 
are  obtained  for  n<m. 


the  trivariate  interferogram.  The  actual  data  fusion 
process,  which  generates  a  single  signal  that  encodes  in  a 
correlated  way  the  three  types  of  data,  consists  of 
sampling  the  trivariate  interferogram  along  a  particular 
sampling  path.  The  characteristics  of  the  fusion  process,  i. 
e.,  the  way  in  which  the  contributions  of  certain  spectral 
details  are  enhanced  relative  to  others,  are  controlled 
through  the  geometry  of  the  sampling  path  and  thus, 
eventually,  through  the  adjustable  instrument  parameters. 


Thus,  the  geometric  interpretation  of  the  direct  transform 
of  trivariate  spectrometry  becomes  the  basis  for 
programming  a  multivariate  M-PEMS  as  a  data  fusion 
optical  processor. 

Furthermore,  the  data  fusion  process  performed  by 
the  M-PEMS  can  be  made  to  adapt  to  the  actual  data  so  as 
to  maximize  its  discrimination  and  identification  power. 
This  can  be  achieved  by  controlling  the  instrument 
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parameters,  and  thus  the  sampling  path,  through  a  closed 
loop  driven  by  an  external  data  postprocessor. 

4.2  External  Data  Postprocessing  for  Discrimination 
and  Identification 

With  suitable  data  acquisition  and  processing 
electronics,  PEM/FT  spectrometers  can  acquire  and 
linearly  extract  spectral  parameters  during  the  very 
scanning  of  the  interferogram,  with  additional  dead  times 
of  less  than  1  ps.  This  is  made  possible  by  the  use  of 
dedicated  real-time  linear  parallel  processors  for  spectral 
deconvolution,  the  first  of  which  was  developed  by  us  in 
the  mid-80s  (Buican,  1985). 

A  linear  processing  stage  can  greatly  reduce  data 
dimensionality  at  the  full  speed  of  data  acquisition.  A 
second,  nonlinear  stage  in  a  data  processing  pipeline  can 
perform  high-performance  discrimination  and 
identification  tasks.  Our  work  on  supervised  and 
unsupervised  neural  network  learning  may  provide  an 
integrated  neural  network  architecture  for  spectral 
identification  (Carrieri,  1999)  offering  distributed  analytic 
capabilities  to  characterize  known  threats  by  type  and 
concentration.  Such  backward-error  propagation  and  self¬ 
organization  mapping  analytic  neural  network  engines 
will  incorporate  technology  capable  of  learning  to 
characterize  IEDs  and  unknown  CB  substances  in  the 
specific  M-PEMS  instrumentation,  data  acquisition, 
preprocessing,  and  software  context. 

5.  PRELIMINARY  DESIGN  CONSIDERATIONS 

M-PEMS  requires  accurate  and  stable  control  of  the 
geometry  of  the  sampling  path  and  thus  of  the 
interferometers’  retardation  amplitudes  and  oscillation 
frequencies,  as  well  as  of  their  relative  phase.  While 
commercial  PEMs  are  still  essentially  free  running 
oscillators  whose  frequency  drifts  uncontrollably  with 
temperature  and  amplitude,  our  work  on  driving  and 
controlling  PEMs  and  PEM  stacks  at  high  retardation 
amplitudes  (Buican  and  Carrieri,  2004;  Buican,  2005, 
2006)  has  eliminated  this  problem. 

We  showed  (Buican,  1997)  that  spectral  resolutions 
of  10  cm'1  and  better  can  be  achieved  in  the  infrared  with 
stacked  PEMs.  The  time  (delay)  resolution  can  be 
estimated  from  Equations  (5)  and  (8)  as  being 
Sp « Scox /(co0xcox)  and  thus  dependent  on  the  excitation 

optical  frequency,  as  well  as  on  the  frequency  of  the 
mechanical  oscillation  of  the  excitation  PEM-based 
interferometer.  For  a  spectral  resolution  of  10  cm"1,  an 
excitation  wavelength  of  1  pm,  and  a  PEM  frequency  of 
30  kHz,  the  resulting  delay  resolution  is  of  the  order  of  5 
ns.  In  a  lidar  application,  this  corresponds  to  a  range 
resolution  of  0.75  m.  It  should  be  noted  that 
superresolution  can  be  achieved  when  the  number  of  basis 


spectra  is  limited  to  those  associated,  for  example,  with 
known  threats. 

We  also  showed  in  subsection  3.3  that  the  duration  of 
a  full  scan  of  the  sampling  path  is  proportional  to  the 
density  of  the  coverage  of  the  sampling  cylinder  by  the 
path.  Preliminary  results  suggest  that  fairly  low  values  of 
the  scanning  parameters  m  and  n  such  as  n/m=6/ 5  result 
in  response  functions  that  are  likely  to  suffice  in  many 
important  applications.  With  an  emission  interferometer 
running  at  30  kHz,  this  results  in  the  acquisition  of  6,000 
double  interferograms  per  second  and  a  trivariate 
spectrum-based  analysis  roughly  every  84  ps. 

A  multivariate  PEM/FT  spectrometer  can  cover  the 
spectral  region  from  the  UV  to  near  IR  with  one  set  of 
optics  (fused  silica),  and  from  the  red  to  12  pm  with 
another  (ZnSe).  Thus,  the  advantages  of  intrinsic 
multivariate  data  fusion  are  complemented  in  PEM/FT 
spectrometers  by  extended  spectral  range,  which  makes 
possible  discrimination  and  identification  based 
simultaneously  on  a  very  wide  range  of  mechanisms  of 
interaction  between  targets  and  optical  radiation. 

6.  APPLICATIONS  FOR  IED  AND  CB  DEFENSE 

Multivariate  PEM/FT  spectrometry  is  potentially  the 
enabling  technology  for  fast  multivariate  biological  and 
chemical  detection  and  mapping.  We  make  this  claim 
because  even  in  its  most  basic  implementation  an  M- 
PEMS  instrument  offers  a  more  complete  description  of 
the  interactions  between  light  and  sample/target  than  any 
other  single  instrument.  Indeed,  in  basic  M-PEMS,  the 
target’s  interactions  with  light  are  represented  by  a 
trivariate  spectrum  in  which  spectral  properties  in  terms 
of  excitation,  emission,  and  delay/lifetime  are  depicted  in 
a  correlated  way.  When  used  for  discrimination  based  on 
tri variate  spectral  signatures,  the  interferogram  signal 
carries  encoded  information  allowing  for  discrimination 
on  the  basis  of  spectral  differences  in  this  3D  space  of 
spectral  variables.  In  other  words,  the  same  M-PEMS 
system  can  not  only  simultaneously  discriminate,  for 
example,  between  different  fluorescent  species  or 
between  elastic  scatterers  at  different  ranges,  but  can  also 
simultaneously  resolve  fluorescence  from  inelastic 
scattering,  while  extracting  some  range  information  as 
well. 

The  basic  M-PEMS  technology  (i)  provides  active 
spectral  analysis,  without  having  to  rely  on  external  and 
uncontrollable  light  sources;  (ii)  is  broadband,  with  only 
two  optical  materials  covering  VUV  to  NIR  and  NIR  to 
12  pm,  and  thus  covering  electronic,  vibrational,  and 
rotational  transitions  and  their  interactions;  and  (iii) 
operates  at  ultra-high  speeds,  with  an  interferogram  being 
scanned,  acquired,  and  processed,  depending  on 
instrument  configuration,  in  an  interval  from  several 
microseconds  to  milliseconds — it  can  thus  deal  with  fast 
changing  targets  in  single-shot  mode,  and  potentially  with 
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ultrafast  light  sources  and  target  events  in  a  repetitive 
excitation  mode. 

To  these  interesting  capabilities  one  can  add  the 
eminently  feasible  concept  of  an  ultra-broadband  M- 
PEMS  instrument  covering  essentially  the  full  spectral 
range  from  the  UV  to  12  pm.  With  suitable  light  sources, 
such  an  instrument  could  simultaneously  analyze  (i)  the 
UV-VIS  LIBS  fluorescence  of  ions  and  neutral  species  in 
the  plasma  plume  created  by  an  ultra-short  laser  pulse;  (ii) 
the  possible  fluorescence  of  molecular  fragments  at  the 
periphery  of  the  ablation  zone;  and  (iii)  the  IR  emission  of 
the  same  molecular  fragments.  Such  a  system  would 
complement  the  abundant  but  potentially  ambiguous 
LIBS  data  with  the  much  more  specific  molecular  IR 
spectral  information.  Such  a  system  would  likely  be 
useful  in  identifying  specific  surface  contaminants 
associated  with,  for  example,  improvised  explosive 
devices  (IED). 

The  extension  of  multivariate  spectrometry  to  include 
the  spectral  dependence  of  changes  in  polarization 
resulting  from  the  interaction  between  target  and  the 
modulated  illumination  beam  promises  to  further  increase 
the  discrimination  and  identification  capabilities  of 
multivariate  spectrometry.  Indeed,  work  done  on 
differential-absorption  Mueller  matrix  spectroscopy 
(DIAMMS)  (Carrieri  et  al.,  1998;  Carrieri,  1999) 
indicates  that  this  may  well  be  the  case. 

If  this  potential  is  realized,  M-PEMS  can  become  the 
technique  of  choice  for  many  important  applications 
aimed  at  protecting  the  warfighter  on  the  battlefield  and 
the  civilian  population  in  the  homeland.  Potential 
applications  of  M-PEMS  range  from  early  detection  and 
reliable  identification  of  biological  agents  and  IEDs,  to 
compact  and  fast  detectors  built  into  smart  munitions  that 
could  recognize  contaminated  targets,  as  well  as  to  fast 
automated  CB  and  IED  defense  systems.  The  civilian 
population  would  also  benefit  from  the  widespread  use  of 
such  highly  sensitive  and  reliable  detectors  of  biological 
agents,  which  could  be  used  to  monitor  in  real  time  large 
public  spaces  such  as  airports  and  supermarkets,  or 
confined  spaces  such  as  HVAC  air  ducts. 

CONCLUSIONS 

We  believe  that  active  multivariate  PEM/FT 
spectrometry  has  the  potential  to  become  the  enabling 
technology  for  fast,  sensitive,  and  reliable  optical 
detection  and  identification  of  specific  chemical  and 
biological  agents  and  contaminants.  We  also  believe  that 
active,  correlated,  multivariate  spectrometry  is  superior  to 
other  optical  techniques  in  that  it  offers  comprehensive 
fused  optical  data  in  a  form  that  is  suitable  for  the  highly 
efficient  and  reliable  detection  and  identification  of 
threats. 

By  providing  a  comprehensive  correlated  description 
of  optical  properties  that  can  only  be  richer  in  information 


than  the  more  common  uncorrelated  descriptions,  M- 
PEMS  may  open  a  new  window  on  the  optical  properties 
of  complex  biological  materials.  Independently  of  any 
security  applications,  M-PEMS  thus  has  the  potential  to 
become  the  spectrometric  technique  of  choice  for  the 
study  of  dynamic  biological  and  chemical  systems. 
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